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Abstract. Distances to AGB stars with optically thick circumstellar shells cannot 
be determined using optical parallaxes. However, for stars with OH 1612 MHz maser 
emission emanating from their circumstellar shells, distances can be determined by the 
phase-lag method. This method combines a linear diameter obtained from a phase-lag 
measurement with an angular diameter obtained from interferometry. The phase-lag of 
the variable emission from the back and front sides of the shells has been determined 
for 20 OH/IR stars in the galactic disk. These measurements are based on a monitoring 
program with the Nan 9 ay radio telescope ongoing for more than 6 years. The interfero¬ 
metric observations are continuing. We estimate that the uncertainties of the distance 
determination will be ~20%. 


1. Introduction 


In the final phases of AGB evolution, the mass-loss rate increase is such that the star is 
partially obscured in the visual and even near infrared bands, as the circumstellar enve¬ 
lope (CSE) gets optically thick. It is assumed that this occurs mostly for intermediate 
mass stars ( M > 2 M 0 ), but direct evidence is lacking, at least for the Galactic disk, 
because of uncertainty in the distances. 

Distances to obscured AGB stars cannot be determined using optical parallaxes; 
however, the maser emission shown by many OH/IR stars can be used. Parallax mea¬ 
surements of pos t-AGB stars and Mira variables were successfully made using H 2 O 
( Tmai et all 120071) and OH masers dVlemmings & van Langevelde! 120071) . but have not 
been attempted for OH/IR stars so far. Another promising way to obtain reliable dis¬ 
tances is the measurement of the phase-lag between the varying peaks of the OH maser 
emission. This yields the absolute linear diameter of the shell and, in combination with 
the angular diameter of the shell measured in interfero metric maps, gives the distance 
( Herman & Habin3 1 1 9851: van Langevelde etahl 1990h . The OH/IR stars usually ex¬ 
hibit maser spectra dominated by two peaks with a separation of 30 ± 10 km s' 1 . This 
double-peak characteristic can be explained if the maser emission comes from a rather 
thin, uniformly expanding shell. The strongest masers then originate from the front and 
the back sides of the shell. The stars are often large-amplitude variable stars, and the 
maser flux density follows these regular variations. Due to the light travel time through 
the CSE, the observer sees the peak from the back side of the shell responding with a 
delay tq relative to the peak on the front side. These delays are of the order of weeks. 
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With the Nan 9 ay Radio Telescope (NRT) we monitored a sample of 20 OH/IR 
stars to measu re their phase-lags. For ab out half of the sample, ph a se-lag s had been 
determined by iHerman & Habind dl985h and Ivan Langevelde et all (Il99dll . Some of 
the stars served as comparison stars to verify the reported phase-lags. For others a 
re-determination was the goal, as the reported phase-lags were inconsistent with each 
other. For the remaining stars, phase-lags were determined for the first time. We have 
also started to image stars at 1612 MHz with the eVLA and e-MERLIN near maximum 
light to determine the angular diameters. We report here on the status of the program 
and give preliminary results. 




Figure 1. Sample light curves of the OH maser flux integrated over the full velo¬ 
city profile for IRAS 20234-1357 (P = 1.16 yr) and OH 30.1-0.7 (P = 5.95 yr). 


2. Light Curves and Phase-lag Determinations 

The sample of OH/IR stars was observed with the NRT in the 1612 MHz OH maser 
line once every month over 5 years (2008-2012). For several sources with insufficient 
coverage of the variation cycle, the monitoring is continuing. Typical integration times 
on source were 5-10 minutes, yielding a typical noise level of 100 nrJy. The spectra 
were baseline fitted and calibrated (nominal error ~ 5%) and have a velocity resolution 
of 0.035-0.070 kms' 1 . Representative light curves of the integrated flux arc shown in 
Figure [Q All stars are large-amplitude variables with periods P ~ 1-6 years. The 
light curves are smooth with some month-to-month scatter. This scatter is not due to 
secular variations of individual maser features in the spectral profile, superposed on the 
smooth variation of the flux integrated over all features. Indeed, dividing spectra of 
adjacent months showed that any variations affected the full profile, and not individual 
features only. We therefore conclude that the month-to-month scatter is dominated by 
calibration errors. 

The light curves of objects with periods < 2 years are well represented by a sine 
curve, except that the level of maximum light is not necessarily recovered in each cycle. 
The light curves become increasingly asymmetric with longer periods, with a steeper 
rise towards maximum (duration 0.3 - 0.4 • P ) and a gentle decline towards minimum. 
Simple analytical functions, such as an asymmetric sine function, are not able to repro¬ 
duce the detailed shapes of these light curves. 
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To find the phase-lag To between the front- and back-side emission (corresponding 
to the blue and red parts of the maser velocity profile), we assumed that the maser 
light curves of both velocity ranges are similar, e.g. they differ only in mean flux and 
amplitude. For all sources, the peak flux densities of the two maser peaks were always 
»1 Jy, so that all parts of the light curves were sampled with high signal-to-noise ratio. 
The phase-lags were then determined by minimizing the function A F = F r {t)-a-Ff,{t + 
to) + c, where /-),. F, are the integrated fluxes over the blue and red parts of the velocity 
profile, and a and c are constants to match amplitude and mean flux of the two maser 
features. An example is shown in Figure O for the OH/IR star OH 16.1-0.3, which 
has one of the longest periods in the sample: 6.03 years. In this case the uncorrected 
light curves show an apparently smaller phase-lag in the rising part of the light curve 
and a larger lag in the declining part. This case shows that full coverage of a variation 
cycle is essential to determine amplitude and mean flux correctly before the phase-lag 
is determined. The resulting phase-lag of 100 ± 10 days is more than twice as large as 
the value, 38.3 days, measured by Ivan Langevelde et all (tl990h . although their period 
of P = 6.2 years is in good agreement with ours. 




Figure 2. Left: Light curves of the maser emission integrated over the blue part 
([-3;+21] km s , crosses) and the red part ([+21;+45] km s _1 ; circles) of the 
velocity profile of OH 16.1-0.3 (P = 6 yrs). Right: Best match of the light curves 
after scaling and shifting the blue light curve and applying a phase-lag to = 100 
days. 


3. Improved Angular Diameters 

The determination of the angular diameter in principle requires imaging the OH maser 
shell at tangential velocities, e.g. in the center of the velocity interval, where the emis¬ 
sion is faint. Because this emission is hard to detect, the diameters are usually deter¬ 
mined by extrapolation of the projected diameters at non-tangential velocities assuming 
that the masers reside in a thin radially symmetric shell. We are currently observing the 
stars interferometrically, close to the maxima of their light curves, to detect as much 
as possible from the faint tangential emission. The aim is to improve the angular dia¬ 
meters derived from less sensitive observations in the past, to determine diameters for 
stars not imaged at 1612 MHz so far, and to derive constraints on deviations from radial 
symmetry of the shells. 
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OH 16.1 -0.3 is one of the objects we observed with the eVLA, and which had not 
previously been imaged at 1612 MHz. The interferometric observations were obtained 
in November 2011 with a velocity resolution of 0.69 km s _1 . Calibration, following 
the relevant steps for observations taken in spectral-line mode, and final imaging were 
performed with CASA (Common Astronomy Software Application). A sensitivity of 
2 mJy beam -1 per channel was reached. The restoring beam size was (T/95 x T.'O) with a 
position angle PA « -28°. Figure [^presents, in the right panel, the integrated emission 
over the red inner wings of the spectra covering the velocity range [+25; +37] km s _1 , 
which gives the best rendering of the diameter of the shell, and a typical spectrum of 
the source at the left. The shell is roughly circular with a possible hint of an elongated 
streak of material in the south east of the shell. 

The best fit (by eye) of the projected shell diameter is given by the magenta circle 
of ~ 3743 diameter. Correcting for the projection effect assuming a thin shell, the full 
angular diameter of the shell is tf> = 471 ± 076. In combination with the linear diameter 
r/oH = (1-73 + 0.17) 10 4 AU corresponding to the phase-lag of 100 days, the distance to 
OH 16.1-0.3 is D = 4.2 ± 0.7 kpc. 



Figure 3. Left: NRT spectrum of OH 16.1-0.3. Right: Integrated emission over 
the red inner part of the spectrum of OH 16.1-0.3 covering the velocity range 
[+27; +37] km s _1 . The magenta circle presents the best fit of the projected diameter 
of the shell. 


4. Results and Conclusions 


Preliminary results of the observing program of the full sample are summarized 
in Table [T] Given are the object name, the newly determined period P, the linear 
diameter don (in AU) = 0.1731 ■ to (in days) calculated from the phase-lag to, the 
angular diameter and the distance D. Uncertainties are < 3% for the periods P 
and approximately 20% for the linear diameters don- P and don are taken from the 
NRT observations, while the angular diameters of the shell are taken from the literature 
( Amiri etall 20121; Baudlll98ll:1chapman et al.lfl984l; iHerman et al.ll 19851: iNorris et al. 


1982; Wolak et al. 


20131) . except for OH 16.1-0.3 (eVLA) and for OH 83.4-0.9 (e 


MERLIN), which come from our own recent observations. 

The shell diameters cover a range of ~ 2000 - 17 000 AU, and in general the dia¬ 
meters are larger for longer periods. The angular diameters taken from the literature can 
be rather uncertain, depending on the degree of extrapolation applied to derive the full 
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extent of the shells. An obvious case is OH 141.7+3.5, where the angular diameter pro¬ 
vided by IChapman et al.1 ( 19841) is strongly underestimated, yielding an unreasonably 
great distance. The accuracy of the distances currently is not better than 20%, where 
this uncertainty reflects only the uncertainty of the phase-lag determination. 


Table 1. Results from the OH 1612 MHz monitoring program. 


Object 

P 

[yd 

rfoH 

[10 3 AU] 

<P 

["] 

D 

[kpc] 

Object 

P 

[yr] 

doH 

flO 3 AU] 

</> 

["] 

D 

[kpc] 

20234-1357 

1.16 

< 1.7 



OH 32.0-0.5 

4.16 

11.2 

1.07 

10.6 

OH 44.8-2.3 

1.47 

2.6 

2.6 

1.0 

OH 127.8+0.0 

4.36 

9.5 

2.69 

3.6 

IRC +50137 

1.74 

1.7 



OH 26.5+0.6 

4.36 

6.0 

4.49 

1.4 

01037+1219 

1.78 

4.3 

8.0 

0.5 

OH 75.3-1.8 

4.52 

10.4 



05131+4530 

2.88 

5.2 



OH 32.8-0.3 

4.63 

10.4 

2.74 

3.8 

OH 39.7+1.5 

3.45 

2.9 

4.0 

0.7 

OH 20.7+0.1 

4.71 

15.6 

1.65 

9.4 

OH 55.0+0.7 

3.48 

8.6 



OH 104.9+2.4 

4.79 

6.0 

2.88 

2.1 

21554+6204 

3.51 

5.2 



OH 30.1-0.7 

5.95 

10.4 

2.69 

3.9 

OH 138.0+7.2 

3.86 

1.7 

0.8 

2.2 

OH 16.1-0.3 

6.03 

17.3 

4.10 

4.2 

OH 83.4-0.9 

4.11 

5.2 

~1.8 

-3.0 

OH 141.7+3.5 

6.05 

13.0 

0.76 

17.1: 


For eleven stars of the sample, phase-lags (linear shell diameters) are available 
from van Langevelde et al. dl990l) . Within our estimated uncertainty of 20% we ob¬ 
tained agreement for OH 26.5+0.6, OH 32.0-0.5, OH 44.8-2.3 and OH 104.9+2.4. 
For the remaining sources, as in the case of OH 16.1-0.3, higher deviations were found. 
OH 127.8+0.0 was monitored bv lWolak et al.l(i201 3l) between 2002 and 2008, also with 
the NRT. They determined a period of 1600 days and a phase-lag of 62 days, which is 
in perfect agreement with our results (P = 1590 days and to = 55 days). 

We expect that the linear shell diameters can be determined with an error ~ 10%, 


by analyzing the variations of the phase-lag across the velocity profile (see Wdlak et al. 


2013]), instead of using the integrated fluxes. In addition, the uncertainties in the dis¬ 
tances will depend on the deviations from radial sym metry of the shells. Assuming 
that these deviations do not exceed 10% dEtoka & Diamond! l2010l) . the uncertainties 
in the distances should not exceed 20%. This will be probed by the interferometric 
observations, which are expected to continue until 2017. 
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Discussion 

Whitelock: Very nice results. Do you have accurate luminosities from the mi d-infrared 
to go with the distances? 

Engels: Bolometric fluxes can be obtained nowadays between 2 and 200 gm by various 
IR surveys. Combining them with the phase-lag distances gives a range of 5000 - 
50000 Lq, consistent with a classification as intermediate-mass stars at the very end of 
AGB evolution. 

Olofsson: If your stars arc in the super-wind phase, it appeal's to me that you find much 
larger durations of this phase than the ones obtained from dust studies (see Lombaert’s 
presentation). 

Engels: The masers are pumped by photons emitted in the superwind at much smaller 
distances from the stars than the masers themselves. The masers might be located at 
distances where the dust and gas is from the pre-superwind phase. 



